Electron-phonon interaction in the solid form of the smallest fullerene C20 
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The electron-phonon couphng of a theoreticaUy devised carbon phase made by assembhng the 
smallest fuUerenes C20 is calculated from first principles. The structure consists of C20 cages in an 
fee lattice interlinked by two bridging carbon atoms in the interstitial tetrahedral sites (/cc-C22)- 
The crystal is insulating but can be made metallic by doping with interstitial alkali atoms. In the 
compound NaC22 the calculated coupling constant A/A''(0) is 0.28 eV, a value much larger than in 
^ . Ceo, as expected from the larger curvature of C2o- On the basis of the McMillan's formula, the 

calculated A=1.12 and a /i* assumed in the range 0.3-0.1 a superconducting Tc in the range 15-55 
K is predicted. 
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PACS numbers: 71.15.Mb, 74.70.Wz, 74.10.- 



O \ The discovery of high superconducting transition temperature in alkali-doped Cgo (Tc= 10-40 K) has attracted much 
I attention on these materials as representatives of a new class of superconductors [Q, ^ . Superconductivity, albeit at 
' lower temperature (Tc=7 K), has been recently observed in another member of this class, C70, which indicates the 
possibility that other fuUerenes might also become superconducting j^. In this respect, the smaller fullerenes are the 
most interesting since the electron-phonon coupling (and hopefully Tc) has been calculated to become larger with 
•1^ . increasing curvature of the molecule |[ ^, |^. Therefore, the recent synthesis of the smallest fullerene, having the 
^ ' largest curvature, C20, starting from the stable C20H20 molecules Q represents a step further along the way towards 
^ high-Tc superconductivity in carbon-based materials. Experimental evidence for the realization of the next step, the 
formation of a crystalline solid based on C20 clusters, has been recently provided as well Crystalline micro-domains 
of the new carbon phase have been detected by Trasmission Electron Microscopy (TEM) in carbon films deposited 
Q on a nickel substrate by the ultraviolet laser ablation of diamond. A face-centered-cubic {fee) symmetry have been 
Q ^ inferred from TEM data. The Raman spectrum shows a strong peak at 1618 cm^^ typical of the stretching of the 
. ethylenic C=C double bond and the presence of C20 clusters as building blocks of the new crystal has been determined 
' by mass spectroscopy of laser-desorbed fragments Based on this information we have theoretically conceived a 
^ . C2o-based solid whose calculated structural and vibrational properties are suitable to reproduce the experimental 
TEM and Raman data . The structure comprises C20 cages in a fee lattice interconnected by two bridging atoms 
00 per unit cell in the tetrahedral interstitial sites (Fig. |l|). Hereafter it is referred to as /CC-C22, the unit cell containing 
^""^ ' 22 atoms. 

[ In this letter, we report on the ab-initio study of the electronic properties of /CC-C22 doped with alkali metals (Na, 
, Li) in the search for metallization and superconductivity. The calculated electron-phonon coupling constant A for the 
' alkali-doped /CC-C22 is found to be substantially enhanced compared to Cgo leading to the possibility of larger than 
-j~>; in alkali-doped fuUerite. 

We use density functional theory in the local density approximation (LDA) , norm-conserving pseudopotentials and 



^ \ plane- wave expansion of the Kohn-Sham orbitals up to a kinetic cutoff of 40 Ry |11 . Non-linear-core corrections 
I have been added for the alkali ions. Phonon frequencies and electron-phonon coupling constants have been computed 
' O within the linear response theory ||l^, |l^. This approach has been successfully applied to the calculation of phonon 
^ . dispersion relations in insulator and metals p2| , |l4j and of the electron-phonon coupling constant A in simple and 
^ ■ transition metals |15 . Brillouin Zone (BZ) integrations are performed using Hermite-Gaussian smearing of order one 
• • with a Hnewidth of 0.01 Ry and 18 special points in the irreducible wedge of the BZ (IBZ). 

^ ■ The calculation of the equation of state of /CC-C22 yields a cohesive energy of 0.63 eV/atom lower than that of 
diamond and a lattice constant of 8.61 A. The density is 1.78 g/cm^, close to that of fuUerite (1.68 g/cm'^), but 
5J] ■ the bulk modulus is as large as 239 GPa due to the presence of covalent interlinks between adjacent C20 cages. The 
' crystal belongs to the Fm3 space group and contains only three independent atoms in the unit cells, the other being 
obtained by symmetry operations. Their coordinates in units of the lattice constant are A = (0.25, 0.25, 0.25), B = 
(0.148, -0.148, 0.148), C = (0.078, -0.217, 0.0) 0, [l6| (Fig. 1). Eight atoms out of twenty of the cage (atoms B) 
and the interstitial carbon atoms {A) are sp'^ hydridized. The other twelve atoms of the cage (C) are sp^ hybridized 
forming C=C ethylenic bonds. There are three bond lengths in the crystal: the single sp'^-like B — C bond in the 
pentagonal ring (1.531 A), the C=C ethylenic-like bond (1.346 A) long and the single sp'^-like bond between a cage 
(equivalent to B) and an interstitial atom {A) (1.517 A). 

The electronic band structure of /CC-C22 is reported in Fig. ^j. The system is an insulator with a slightly indirect 
XL band gap of 2.47 eV. Projection of the Kohn-Sham states on the minimal atomic basis set shows that the lowest 
conduction bands and the highest valence bands are mainly formed by tt* and tt states of the ethylenic-like bonds. 
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It is suggestive to wonder if doping might turn /CC-C22 metallic by electron transfer into the molecular-like lowest 
conduction bands, as occurs for instance in fuUcrite. The analogy with Cgo is evident by considering that both Cgo 
molecules in fuUerite and C20 cages in /CC-C22 form a fee lattice. In the A3 Cgo (A=K, Cs, Rb) compounds the 
three alkali metal ions per unit cell occupy the two tetrahedral sites and the octahedral site 0. In /CC-C22 the two 
tetrahedral sites are occupied by the interstitial carbon atoms, but the octahedral site is empty and might host a 
guest ion. In the search for metallization we have inserted a Na atom per unit cell in the octahedral site. Calculation 
of the equation of state of the compound NaC22 shows an increase in the lattice parameter (a=8.765 A) with respect 
to the pure phase and a larger (2 %) increase in the C=C bond length due to population of antibonding tt* states. In 
fact, Na is ionized by donating its outermost electron to the host which becomes metallic. The Na-host interaction 
is weak and the band structure of NaC22 is very close to that of the pure system, but for the shift of the Fermi level 
inside the lowest conduction bands (Fig. |^b). Therefore the rigid band approximation holds to a large extent in the 
description of Na insertion in fcc-C22- Although the lowest conduction band is 3-fold degenerate at the F point (Fig. 
Ija) , the Fm3 symmetry is preserved in NaC22 and no static Jahn- Teller distorsion is observed. This is due to the fact 
that the minima of the lowest conduction bands are located at the X and L points where the degeneracy is already 
lifted. 

To study the superconducting properties we have computed the electron phonon coupling constant A as 
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where iV(0) is the electronic density of states at the Fermi level (zero of energy) per spin and per cell, Aaq/iV(0) 
is the partial electron-phonon coupling due the a-th phonon at the point q in the BZ (of volume f^Bz), Caq is the 
phonon polarization vector, M is the atomic mass matrix, n and m are electronic band indices, Uk.n is the periodic 
part of the Kohn-Sham state with energy i?kn, and VVjjy is the derivative of the Kohn-Sham effective potential with 
respect to the atomic displacement caused by a phonon with wavevector q. We replaced the (5-functions in Eq. |^ 
with Hermite-Gaussian forms of 0.01 Ry linewidth. The BZ integration over k is performed on a grid of 18 special 
points in the IBZ. We computed Aaq/7V(0) for q along the high symmetry directions [100] and [111]. The integration 
over q in Eq. ^ is performed over a uniform mesh in the IBZ by assigning to each point the value of Xaq/N{0) 
computed for the closest point along a high symmetry directions [0. Finally, for NaC22 we obtain A/A^(0)=0.28 eV, 
N{0)—4:.0 states/eV/spin/cell (as computed by the tetrahedron method with 213 points in the IBZ) and A=1.12. A 
similar calculation for the charged fuUerite Cqq at the F point only in the molecular-like approximation gives 
A/A^(0)=0.076 eV, in agreement with previous works [^ (or 0.064 eV by including the set of phonons considered in 
table IV of Gunnarsson's review [Q). Therefore, although /CC-C22 is not strictly a fuUerite, the C20 clusters being 
covalently interlinked via interstitial atoms, the solid retains some features of its molecular building blocks, giving rise 
to a substantial enhancement of A/A^(0) with respect to Cgo as expected from the larger curvature of C20 [^. The 
persistence of molecular features in the crystal is also confirmed by the calculation of A/A^(0) in the molecular- like 
approximation [|l^ for the pure /CC-C22 doped with an additional electron (equally shared by the three lowest tu 
states). This latter calculation gives A/A^(0)=0.31 eV, a value very close to the result for NaC22 at convergence in BZ 
integration (0.28 eV) |l^. A tight-binding calculation of the electron-phonon interaction in the isolated C20 cluster 
gave a somehow lower value A/A^(0) = 0.18 eV Q which may be partially due to the semiempirical method used. 

In order to identify the phonons which mostly contribute to the elctron-phonon coupling, we report in Table I 
the partial electron-phonon coupling constants Aq/A'^(0) of the individual a-th phononic band, obtained by averaging 
Xaq/N(0) over the BZ. The frequency and symmetry of the a-th phononic bands at the F point are also given in Table 
I for the pure /CC-C22 and NaC22 phases. Although the phonons with the largest matrix elements in Eq. ^ correspond 
to the stretching modes of the C=C ethylenic bonds (1359-1362 cm~^), their contribution to A is attenuated by the 
factor l/cj^. As a consequence the largest Xa/N{0) is due to the Eg mode at 334 cm~^ (cfr. Table I); its displacement 
pattern (Fig. ||) produces the largest modulation of the dihedral angles which control the tt conjugation of the ethylenic 
bonds. This is at variance with Cgo where the largest electron-phonon coupling constants are due to the phonons 
with the highest frequencies 

In system where the superconductivity is driven by the electron-phonon interaction an estimate of Tc can be obtained 
by making use of the McMillan's solution [^0[ |l| of the Eliashberg equations: 

= T^exp , (2) 



1.2 A-^*(1 + 0.62A) 
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where ujin is the weighted logarithmic average of the phonon frequencies [^0[ (470 cm~^ in NaC22) and /i* is the 
Coulomb pseudopotential n = UN{Q) renormalized by retardation effects. J7 is a typical screened Coulomb interaction 
among states at the Fermi level |2^. Latest estimates give fj,* ^ 0.3 for Cgo although, due to the similarity between 
the electronic bandwidth (0.5 eV) and phonon energies (up to 0.2 eV), the applicability of Eq. || is questionable 
for fuUerite j^] . We have not attempted to estimate numerically fi* for NaC22, so our discussion now turns on a 
very speculative level. Although U might be larger in NaC22 than in fuUerite due to the smaller size of C20, ^(0) 
turns out to be lower (4.0 and 6.6-9.8 states/eV/spin/cell in NaC22 and K3C60 from band structure calculations [0, 
respectively) due to the larger dispersion of the conduction bands, three times larger than in Cgo. The estimate of 
the renormalization of /x leading to /i* is difficult , however the renormalization is perhaps larger than in Ceo due 
to the larger bandwidth. Thus, if we assume that the value of 0.3 estimated for fuUerite is an upper bound for ^* in 
NaC22 we obtain Tc=15K from Eq. ^ and A=1.12. For the more optimistic case of fj,* in the range 0.1-0.2 typical of 
simple and transition metals ||2^, |2^, we obtain the remarkable result of Tc=55-33 K {Tc changes linearly with /i* in 
this range). 

Since A/iV(0) is mainly a molecular property, we have also attempted to increase A by increasing iV(0) with a lattice 
expansion. This has been achieved by substituting C by Si in the interstitial sites producing the compound NaSi2C2o. 
This latter system has a substantially expanded lattice parameter (a=9.49 A, A^(0)=4.44 states/eV/spin/cell) but 
only a marginally enhanced A=1.49 with respect to NaC22. In fact the increase in N{0) is lower than expected from 
the lattice expansion due to the reduction of bands at the Fermi level from two in NaC22 (cfr. Fig. H) to one in 
NaSi2C2o. Moreover, a decrease in ujin due to the lattice expansion compensates the increase in A and finally yields a 
Tc close (within 10 K) to that of NaC2o. However, further increase in A is foreseeable for the compounds NaGe2C2o 
and NaSn2C2o. 

In summary, we have carried out ab-initio calculations on the electronic and superconductive properties of a new 
carbon phase /CC-C22 formed by assembling C20 clusters in a fee lattice. The system is insulating but can be made 
metallic by doping with interstitial alkali metals. Although the structure is not a molecular crystal, the C20 cages 
being covalently interlinked via interstitial carbon atoms, the crystal retains some molecular features. In fact, as 
expected from the molecular properties of the isolated C20 cluster, the compound NaC22 exhibits a large electron- 
phonon coupling constant A/A^(0) = 0.28 eV, strongly enhanced with respect to Cgo due to the larger curvature of C20. 
Estimate of N{0) from band structure calculation yields A =1.12 which allows a rough estimate of Tc from McMillan's 
formula. For fj,* in the reasonable range of 0.3-0.1 we predict Tc=15-55 K for NaC22. Although a more precise estimate 
of Tc requires the difficult calculation of fj,* , the results presented here shows that NaC22 is a potentially interesting 
superconductor. We hope that these theoretical results could stimulate further experimental work on the synthesis 
and characterization of this new structure. 
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TABLE I: The partial electron-phonon coupling constant \a/N{0) (eV) of the individual a-th phononic band, obtained by 
averaging Aaq/A'^(0) over the BZ. The frequency and simmetry of the a-th phononic bands at the F point are also given for the 
pure /CC-C22 and NaC22 phases. Only the phononic bands with Aq/A'^(0) larger than 5 meV are reported. Tu(l) correspond to 
the acoustic bands. Tu(2) correspond to translational modes of Na ions. 



Modes Energy (cm"^) Aa/A''(0) 
C22 NaC22 



T„(l) 


- 


- 


0.010 


T„(2) 


- 


271 


0.022 


T«(l) 


281 


282 


0.046 


T9(2) 


490 


445 


0.038 


Efl(l) 


492 


334 


0.072 


T„(3) 


576 


482 


0.019 


Tx.(4) 


626 


609 


0.008 


T8(3) 


701 


637 


0.008 


T„(5) 


731 


697 


0.007 


E„(l) 


754 


719 


0.005 


T.(4) 


857 


823 


0.005 


T„(6) 


948 


927 


0.006 


T„(10) 


1559 


1362 


0.010 


E9(3) 


1574 


1359 


0.005 



FIG. 1: The face-centered-cubic structure of the /CC-C22 crystal. The conventional unit cell with four formula units is shown. 
The crystal belongs to the Fm3 space group and the coordinates of the three independent atoms in units of the lattice constant 
are A = (0.25, 0.25, 0.25), B = (0.148, -0.148, 0.148), C = (0.078, -0.217, 0.0). The atoms C and C" and the other 10 of 
the cage equivalent by symmetry are linked by an ethylenic-like bond 1.346 A long. In the NaC22 compound the alkali atom 
occupy the octahedral site of the fee lattice. 
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FIG. 2: Band structure around the Fermi level of the pure /CC-C22 (a) and NaC22 (b) crystals. The zero of energy is the Fermi 
level of NaC22. The top of the valence bands of the two compounds is aligned. The two lowest conduction bands at the F 
points have tu simmetry. The electronic density of state of NaC22 is also reported (c), as computed by the tetrahedron method 
with 213 points in the IBZ. 




FIG. 3: Displacement pattern of the £^9(1) phonon at 334 cm ^ (cfr. Table I). 



